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The 6 term that breaks the Strong CP symmetry is introduced in the two flavors of dynamical 
QCD simulation. 9 is analytically continued to a pure imaginary number to make the probabil- 
ity of Monte Carlo positive. The Neutron's Electric Dipole Moment (NEDM) is measured on 
the ensemble under a uniform and week electric field. Other applications of 6 terms are also 
discussed. 
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1. Introductions 

Since the first attempt to get the non-perturbative prediction of the Neutron's Electric Dipole 
Moment (NEDM) with the Strong CP violating 6 term, 

f 

Se = ^ 22n 2 I ^^P^^F-cpd x = i ^2to P ; (1-1) 

on lattice 18 years ago [jl]], there have been renewed interests on this calculation recently corre- 
sponding to the new experiments and theoretical developments. The NEDM is very sensitive to the 
sea quark mass, as we will see, and the calculation could be much improved thanks to the recent 
various developments for the simulation with dynamical quarks. 

All the calculations aiming for NEDM so far introduce the effect of Sq in terms of reweighting 
technique: the statistical ensembles of QCD vacuum were generated via the Boltzmann weight of 
6=0 case, Prob(£/ ) oc det Nf D[U ]e^ So ^ , and then the effect of Se is treated as a part of observable 
in the measurement stage, (ff) e = Y,u ^~ 5e Prob(t/)/Lj/ le~ 5e Prob(t/). 

In this study, we introduce the CP violation according to S$ in the probability of QCD en- 
semble generation, then measure observables without reweighting technique to see whether the 
new way of calculation has better control over the statistical and the systematic errors. To keep the 
Boltzmann weight positive semi definite, the value of 6 is analytically continued to pure imaginary, 

6 -> -id. (1.2) 

Such analytic continuation was previously explored for CP" -1 models with 6 term [Q, ||] and QCD 
with the chemical potential [0]. 

To measure NEDM, a constant electric field, E, is applied to Nucleon, and the energy, m^, 
dependence on its spin polarization, S, is measured, 

m N (E,S)-m N (E,-S) = 2id N (G)S ■ E + 0(E 3 ) . (1.3) 

We use electric field constituted from the ordinary U(l) gauge field, 

U* M (x) = e iqA ™ {x \ qA™{x) G R , (1.4) 

on Euclidean lattice, oppose to its Minkowski space version used in the original calculation [jl]], 
which is reported to be susceptible to the systematic error induced at the periodic boundary in 
temporal direction [§]. To avoid the boundary effect and also to safely neglect the GiE?) in ( |Q| ) 
at a time, a uniform and weak electric field invented some times ago is implemented. 
By assuming 6 is small, the NEDM , d^id), is expanded in terms of 6, 

d N (e)=d { N l) e + 0(e 3 ) , (i.5) 



within which approximation, the analytic continuation, (1.2), can be trivially reverted, and more 



importantly, the energy splitting (1.3) becomes a real number, which could be easily measured from 
the ratio between polarized propagators of Neutron whose spin are parallel and anti-parallel to the 
direction of the electric field, 

R(t,t src ;E)= lim ^f^ 8 ^ ^ • M 
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In the following sections, the NEDM calculation and preliminary results are mainly detailed. 
We also briefly report another related attempt, restricting the imaginary 6 action to a time slice and 
identifying it as the source term for the flavor-singlet pseudoscalar meson (T]') for disconnected 
quark loop measurements. 



2. Ingredients for NEDM 

Before describing the novel simulation setups, we compare the external electric field method 
employed in this work with another method, the form factor calculation [B, SJ|, in which NEDM is 
calculated from the CP-odd form factor, F^{q 2 ). While the latter method is free from the €?(E 3 ) 



contamination in ( |1.3| ), one needs to insert momentum q to the vector current of the Nucleon's three 
point function, and extrapolate q 2 — > limit. We chose the simpler method, the external electric 
field method, that involves only two point functions with zero momentum, expecting the statistical 
and/or systematic error could be smaller than the other. 

2.1 Dynamical simulation with 6 terms 

All the previous lattice QCD calculations with the strong CP breaking term Sg, (|lTT|), were 
carried out with the reweighting technique as discussed in Sec. [[]. Alternatively, one could put 
Sg into the ensemble probability if Sg becomes real number by 6 being analytically continued 
to pure imaginary, ([Of). There are two options in the implementation of Sg on lattice: one is the 
gluonic definition from Y lX FF(x) with an appropriate smeared/cooled link field so that the resulting 
topological charge, <2 top , become (close to) an integer value. This implantation may be useful 
for a certain application such as precise determination of the Q top distribution in pure Yang-Mills 
theory [JTCj] . Actual simulation coding, however, might likely involve the smeared/cooled gauge 
force term, and results could strongly depend on the details of the smearing/cooling used. Simpler 
implementation would be using the anomalous axial transformation in the continuum theory, 

vr _,g-»»/(2W / ) yj W ^ We -i6Y 5 /(2N f ) (21) 

-► @yf@\jf e Se , (2.2) 

that cancels the Sg from the original action at the cost of CP violating 75 mass term introduced in 
the fermion's Lagrangian for Np = 2, 

&f= L Wi(P + m + im^Y5)Vi +@{Q 2 )- (2.3) 

i=u,d 

We chose the latter implementation because of the absence the ambiguity from cooling/smearing 
By the analytic continuation, ( |1.2| ), the fermion determinant with the 75 mass term in (2.3) remains 
positive: &et(p + m + im^y 5 ) 2 — ► \det(p + m + m^y 5 )\ 2 . 

Our lattice fermion action is the RG improved clover fermion identifying the continuum 
counter parts as follows: 

p + m^D clmer , im-Y5^Km-m c )-Y 5 , (2.4) 

where rrf is the mass point where the pseudoscalar meson becomes massless. 
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We note the relation, (2.2) (2.3), are not be precisely realized for the clover fermion and there 
will be a discretization error which will vanish in the continuum limit. More demanding simulation 
with overlap fermion or domain-wall fermions might be free from such systematic error. 

2.2 Uniform and weak electric field 

The second ingredients for NEDM is the external electric field on the periodic lattice. The 
electric field needs to be weak and constant everywhere including the periodic boundaries as dis- 
cussed in Sec. [l[ Naively descretizing constant electric field of the continuum theory, for exam- 
ple, in z direction, Af M (x) = Et in terms of vector potential, one obtains uf M = exp(iqEt) for 
charge q fermion field. This Abelian link field gives constant electric field, e' F,z = e' qE , every- 
where except at the boundary between t = L t — 1 and t = 0, where the strong opposite electric field, 
exp(—iqE(L t — 1)), is induced. This delta function like field causes a skewed Neutron propagator 
near the boundary By using the constructions in [Q, [J, such systematic error could be avoided 
by adjusting the temporal link, only at the boundary, as 

U™(z,L t -\)=exp(-iqE(L t -\)z) , £ = |^xZ, (2.5) 



where the condition ( |2.5[ ) is required not to induce another opposite field at the corner plaquette, 
including points (z,t) = (L z — l,L t — 1) and (0,0), of the periodic boundaries. We emphases this 
construction has weaker electric field, E ~ O.Olq x Z than the simpler construction with E = 2n/L t 



so that 0(E ) contamination in ( |l.3| ) is less problematic. 



3. Progress report on NEDM calculation 

We now report the preliminary results of NEDM calculation with the setups described in Sec. || 
employed on 16 3 x 32 lattice. The gauge action of our choice is the RG improved gauge action a la 
Iwasaki with j8 = 2.1(a~ 1 ~ 2GeV). Thermalized statistical ensembles of two flavors of dynamical 
clover fermion action (c sw = 1.47) with 6 = 0.4 and 0.0 are accumulated for z = 2,000 in the 
molecular dynamics time unit. The parameter point and the lattice scale were determined using 



the mean field improvement by CP-PACS [11]. The effect of 6 vanishes at massless limit as 



seen in (2.3). As this is a feasibility study, we consciously set the quark mass unphysically large, 



K = 0.1357 or m K /m p ~ 0.85. The HMC code with accelerations a la Sexton- Weingarten, and also 



Hasenbush is used Q12[]. 

By switching on e~ Se = e~ 2t op in the Boltzmann weight, the negative topological charge 
should appear more frequently than the positive charge for 6 > 0. Figure [j] clearly shows such 
reaction to 6. The positive and negative charge equally emerge for 6 = in the Gaussian-like 
distribution, which is shifted in negative direction, (Q) = —7.9, for 6 = 0.4 ensemble. Here we plot 
the distribution of the iff (a 2 ) improved topological charge defined from the gluon configuration, 
Gtop = §<2 lxl — f<2 2xI , where Q lxl and Q lx2 are topological charges calculated from lxl and 



1x2 clover leafs respectively [13]. The lattice is smoothed until Q top becomes saturated close to 
an integer value. Applying 200 steps of APE smearing (cape = 0.45) turns out to be sufficient. 

The polarized Neutron propagator is measured using the code developed in JT4| ] on 100 con- 
figurations under the external electric field, that is a multiple of, E = 2n/(L z L t ) x 3 ~ 0.0368, in z 
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Figure 1: The distributions of g top for 9 = 0.0 (left) and 9 = 0.4 (right). Q top is calculated by ff{a 2 ) 
improved definition on 200 steps of APE smeared (cape = 0.45) lattices. Also the definition only uses lxl 
plaquette is plotted. The distribution for 9 = 0.4 by the reweighting technique is plotted on the left figure to 
show the poor predictability for such sizable 9. 



direction. The quark charges are set as q up = 2/3 and ^down = —1/3. On each lattice, two sources 
for Neutron made of the Gaussian smeared quarks are prepared, and the propagators of Nucleon 
are measured for each sources separately. 



The spin dependent ratio of the propagator, (1.6), is estimated from the propagator averaged 



with the one with the other source points, and also with the one applied the opposite electric field 
with the appropriate combination as follows: 

/ R(t,t src ;E) \V(2A&«) 
Ri (t - t src ) = TT r \ • (3-1) 

Thanks to the temporal translational invariance of the external electric field described in Sub- 



Sec. 2.2, the backward propagator with the proper parity projection, which is the anti-Neutron, 
at t — t src < may also be used in the average to enhance the signal. While a reasonable signal 
for \t — t src \ < 4 can be seen within current statistics, we need more ensemble and measurements 
to check the excited state contamination. We also compare R%{t) using the traditional reweighted 
technique on 6 = lattice with same statistics, and found about a factor of two smaller statistical 
error for the new method. 

Due to lack of the chiral symmetry in clover fermion action and also from 6{Q 2 ) in (2.3). 



6 in ( |2.3| ) may be mixed with different operators. To check the effect of 8 in other quantities 
than NEDM, we calculate the shift of topological charge (Q t0 p)e and pseudoscalar meson mass 
Mps(d) on 6 = 0.4 configuration. For crude estimates of effective 8 value on this ensemble, 
we assume Gaussian distribution for <2top> an d the leading order prediction of chiral perturbation 
theory, Mj s (d) = Mj s (0)cosh(d /N f ) in [(B]]. The ratio of the effective ds obtained by these 



method to the input value, 0.4, turn out to be 0.6(1) from (<2top)e> an d 0.9(5) from pseudoscalar 
mass. 
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Figure 2: The double ratio, R 2 (t - t src ) = Y\t src (left), and it's effective mass log [R(t)/R(t + 1)J (right) 
are plotted. The applied external electric fields are E = +2Eq (red), E = —2Eq (blue), and the appropriate 
average over E = ±2Eq (black). The red and green symbols are calculated using one source points while the 
block symbols is the average of two source results. 



4. Application for 77 ' calculation 

If we generate ensemble replacing Sq by 
Sx = XY,~WysV{x,t s ) , (4.1) 

which is the source term enhancing the negative topological charge density only at a time slice 
t = t s , the resulting configurations could be used to measure quantities involve the disconnected 
quark loop. Here we present a calculation of the flavor-singlet pseudoscalar meson (17') propagator 
as an example. 

The pseudoscalar density at t with the Dirac operator including the A term, Dx is estimated 
stochastically, 

(Tr7 5 D A V,0>A = -^(Try 5 D-^O75^^ , 

(4.2) 

whose left hand side is the 17' correlation function between t s and t omitting <^(A 3 ). The result of 
effective mass on 1,000 configuration of 8 3 x 16, a ~ 0.2fm lattice is shown in Figure || compared 
with the conventional two point function method (black circle). The central value is consistent with 
each other, however, the source method has much rapid error growth by increasing t — t s . 



5. Summary and Conclusion 

We have generated Np = 2 lattice QCD ensembles with the explicit CP violation term, 6 
term. To preserve the positivity of the generating probability, 8 is analytically continued to pure 
imaginary. By applying the uniform and weak electric field on the ensemble, Nucleon's electric 
dipole moment is calculated. While the results are only on two parameter points, the signal to noise 



6 



Dynamical QCD simulation with 9 terms 



Taku Izubuchi 




Figure 3: The effective mass of the flavor singlet pseudoscalar meson as a function of time using the source 
method compared with the conventional two point function method. 



ratio observed is encouraging considering the relatively small statistics. The checks for the excited 
state contamination by increasing statistics, effects of the terms of higher orders in 6 and E, and 
the quark mass dependence are subjects for future study. 

We also explore QCD ensembles with the source for the flavor singlet pseudoscalar density. 
The propagator of the flavor singlet pseudoscalar meson, r\' , is calculated by restricting the imag- 
inary 6 term on a time slice, consistent results with the conventional method are obtained though 
the statistical error turns out to be larger than that from the conventional method. 
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